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Abstract 

Motivated by certain problems of statistical physics we consider a stationary stochastic process in 
which deterministic evolution is interrupted at random times by upward jumps of a fixed size. If 
the evolution consists of linear decay, the sample functions are of the "random sawtooth" type and 
the level dependent persistence exponent 6 can be calculated exactly. We then develop an expansion 
method valid for small curvature of the deterministic curve. The curvature parameter g plays the role 
of the coupling constant of an interacting particle system. The leading order curvature correction to 
9 is proportional to ~ g^^^ ■ The expansion applies in particular to exponential decay in the limit 
of large level, where the curvature correction considerably improves the linear approximation. The 
Langevin equation, with Gaussian white noise, is recovered as a singular limiting case. 
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1 Introduction 



In this work we study the stationary stochastic process (^(i) that obeys the equation 

dm 



dt 



^-A{0+aY.Sit-t,) (1) 



Here a is a positive parameter and the ti are random times distributed independently and uniformly with 
density p; the random term therefore represents white noise, but with a nonzero average equal to pa. 
Hence ^{t) evolves deterministically except for upward jumps of fixed size a occurring at random times. 
We take the systematic "force" A{^) such that it has positive derivative and satisfies ^(— oo) < pa<A{oo), 
which ensures that ^ possesses a stationary distribution. A special case is the linear equation obtained for 
the choice A(^) — /3£_. Our interest is in the first passage time problem associated with a preestablished 
threshold ^ = X. 

More precisely, for some general stationary process ^{t), let Q{T) be the probability that during a time 
interval of length T it stays above a threshold X, given that it was larger than X at the beginning of that 
interval. For many of the common processes in physics Q{T) decays to zero exponentially with an inverse 
relaxation time 9 defined by 

e= lim T"MogQ(T) (2) 

T — >oo 

Both Q{T) and 6 depend on the threshold value X . 

Physicists are interested in the persistence exponents of various stochastic processes because of their 
connection to critical phenomena; there, after an appropriate rescaling of variables, 9 appears as the 
exponent of a power law and is called the persistence exponent. The theory of critical phenomena has 
brought to light the importance of exponents for the classification of physical systems. This has spurred 
theoretical physicists in recent years to attempt to calculate persistence exponents associated with several 
prominent problems of that discipline. The exponents are in each case nontrivial and unrelated to the 
static and dynamic critical exponents of the same problem. 

Many authors have studied processes of zero average and symmetric under sign change of ^. The 
quantity of primary interest is then "the" persistence exponent associated with the threshold X — Q. For 
nonzero X one also speaks of the level exponents. A review article of earlier work, mainly mathematical, 
is due to Blake and Lindsay [Q. Majumdar ||] and Godreche ^ have provided useful reviews of recent 
work, mainly by physicists. Almost all of this work deals with processes ^(f) that are Gaussian. Among 
these the Gaussian Markovian case is easiest to treat. Majumdar and Sire followed by Oerding et al. 
Q and Sire et al. have designed a perturbative method for processes that are Gaussian and close to 
Markovian. Majumdar and Bray Q have set up an e expansion for smooth Gaussian processes in spatial 
dimension d = A — e. Nontrivial persistence exponents have also been identified for such familiar func- 
tions as the solution of the diffusion equation with random initial condition |^ ^. All these persistence 
exponents appear to be extremely difficult to calculate analytically. 

In physical systems the addition of the effects of many degrees of freedom very often leads to Gaussian 
processes. Nevertheless, certain non-Gaussian processes also arise naturally. In this work we study the 
level exponents for the strongly non-Gaussian case of Eq. (]l|) . One example of how closely related processes 
enter physics is through a question associated with the one-dimensional random walk. Let ^(t) be the 
number of steps needed before the walk has visited t distinct sites. Then ^(t) = e~* is a stationary 

process consisting of exponential decay interrupted by upward jumps. The only difference with Eq. (|^) is 
that it leads to a probability distribution of the jump sizes, whereas in (|l|) we take a a fixed parameter. 
Numerical evidence shows that this problem possesses well-defined persistence exponents, but, as in many 
other cases, no way exists to find them analytically. 

Other instances are furnished by statistical physical models that exhibit cluster growth, such as Ran- 
dom Sequential Absorption and percolation theory; if ^(t) is the suitably scaled size of a particular cluster, 
then jumps are due to coalescence with other clusters. As a specific example, let the bonds of a lattice 
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be rendered percolating in a sequential manner pJ| and let i be the instantaneous fraction of percolating 
bonds; define then ^(t) as the size of the cluster connected to the origin. In spatial dimension one it is 
easily shown that an appropriate scaling (which is such that t--s-ooast^ 1) yields again a stationary 
process ^(t) with a probablity law for the jump sizes. 

In order to study the persistence exponents associated with Eq. (|l|) we exploit the following idea. The 
persistence probability Q{T) is determined by the subclass of^{t) that do not cross below X ioi < t <T. 
When the force A{^) is strongly positive, then a ^{t) in the contributing subclass is unlikely ever to rise 
very high above the threshold X. We conjecture, therefore, that we will obtain a good description of this 
subclass by expanding A{^) around the threshold value ^ = X. We convert this idea into an expansion 
procedure. Roughly speaking, the zeroth and the first order of the expansion are determined by A{X) 
and A'{X)/A{X), respectively, i.e., by the slope and the curvature/(slope)^ ratio of the deterministic 
evolution curve. The precise mathematics is slightly more subtle and shows that instead two parameters 
appear, called r and g, whose definition is more complicated. Our theory then yields the level exponents 
9 in terms of r and g in the small g limit. We shall refer to g as the curvature parameter. 

In Section || we write Q{T) as a path integral on all contributing ^{t). The expression resembles 
the partition function of a system of interacting particles in a one-dimensional volume T, with the jump 
times ti,t2, ■ . ■ in the role of the particle positions. We rearrange the path integral in such a way that 
a "noninteracting" contribution appears, characterized by a parameter r, and a remainder due to an 
"interaction potential" V which is a functional of the jump times. Our use of the term "noninteracting" 
does not mean that the V = problem is trivial - it is not -, but merely that it is purely combinatorial. 
We are led to define the parameter r of the noninteracting theory by 



P 



X+a 



X 



m 



pa (^_aA^^^ 



A{X) \ 2A{X) 

This equation shows that r involves not only A{X) but also the full series of its derivatives. 

In Section || we consider the zeroth order, V = 0. It amounts to replacing A(£_) in Eq. (||) by the 
constant pa/r, so that as a consequence ^(i) is piecewise linear with slope —pa/r. All samples of this 
zeroth order process are therefore "random sawtooth" functions. In this order we shall write the level 
exponent as 9o{X). We find 

6»o(r) =pfilogi-i + l) 0<r<l (4) 
\r r r ' 

For r — > 1 the persistence exponent goes to zero; the interpretation of this unphysical effect is that for 
r > 1 the linearization creates a finite probability for ^(i) to escape to +00. 

In Section ^ we consider the interacting theory, 7^ 0. The potential V is determined by A in a way 
described in that section. We are unable to deal with the general case. Instead, we expand y in a series 
of which we retain only the first term, whose coefficient g plays the role of an interaction constant. The 
expression for g is 

A{X + a) - A(X) 

We show that there are at least two limits in which the higher order terms in the series for V are negligible, 
and in which the remaining problem, with only two parameters r and g, can be solved. 
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The Subsections 4.2 and 4.3 are common to both Umits. The Laplace transform K{Cl) of the path 
integral for Q(T) appears to satisfy a recursion relation whose solution is expressed in Eq. as the ratio 
of two infinite series. Explicit evaluation of these series turns out to be a rather formidable task. The 
soluble limits are the following. 

Limit (i). The limit g at r fixed. Eqs. and 1^ show that this corresponds to A'{X) ^ at 



fixed A{X). In Subsection iA we calculate the exponent 0{r,g) in a small g expansion, with the result 



that a nonanalytic correction term to Eq. appears, 

1 / 1 \ 2/3 /Ott \ 2/3 

0(r,.g) =00 (r)+p -(log -j (—5) g ^ 0, r fixed (6) 

Limit (ii). The limit 17 ^ 0, r — > with fixed ratio g/r. Eqs. (||) and (||) show that this corresponds 
to A{X) ^ 00 at fixed A'{X). In this limit the expansions in Eqs. (||) and (|^) may be replaced by their 
first term, which we shall denote by an index 0, 

pa _ aA'(X) 



This limit, considered in Subsection [4.5| , requires separate analysis; nevertheless, the result for 0{r,g) is 
what one also obtains by naively substituting r — and g — go in Eq. (^) . 

The example of greatest interest is the linear equation that prevails for the choice A{^) — (3^. When 
the threshold X becomes large we have tq = pa/ PX and go = a/X. Upon expressing for this case 6* as a 
function of X we arrive at the explicit asymptotic expansion 



In Section [s] we compare analytical results for both limit cases to Monte Carlo simulations of Eq. (Q) . 
Excellent agreement is found. In particular, there is strong numerical indication that the higher order 
terms in the asymptotic expansion go to zero as X — > 00. 

The name Langevin equation is traditionally reserved for equations of type where the random term 
represents Gaussian white noise. In Section we observe that the white noise of Eq. (|l|) becomes Gaussian 
in the limit p —> 00 and a ^ at fixed pa , and that, correspondingly (and after appropriate rescaling 
of variables) Eq. iMj becomes a Langevin equation. Hence our work enables us to pass continuously from 



strongly non-Gaussian to Gaussian noise. In Subsection 6.1 we place ourselves directly in the Gaussian 
limit and determine, via the associated Fokker-Planck equation, the Gaussian persistence exponent 6g for 
asymptotically high threshold; our method is close to the one of Krapivsky and Redner ||l^. In Subsection 
|6.2| we then investigate how the Gaussian limit emerges from the more general approach of Sections 
Section ^ contains our conclusions. 

2 Phase space integral 
2.1 Solution ^(t) 

The solution of Eq. (|l|) is piecewise continuous. In the time interval between two jumps ^{t) evolves 
deterministically according to 

Cit) = fit-ui), te^i<t<ti (9) 

where Ui determines a shift along the time axis and the function f{t), if we choose it such that /(O) = X, 
is obtained from A{£^) by 
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Hence ui acquires the meaning of the ultimate instant of time at which the £th jump should take place 
if £,{t) is to stay above the threshold. The fact that there is a jump of size a on the border between two 
successive time intervals leads to the identity 

f{ti-ut+i)- f{ti-ue)=a (11) 

We shall be more specific now and consider the solution (^it) of Eq. (Q) with initial value ^(0) = ^o- It is 
uniquely specified by the set of jump times < ii < i2 < • ■ ■ • Eq. ([TTI), which is here valid for £ = 1, 2, ... , 
allows one to express u^+i in terms of t( and w^, and, upon iterating, as a function of ii, . . . , and ui. 
Finally, ui may be eliminated in favor of the initial value by means of /(— ui) = ^o- Hence we have 
obtained the formal answer to the question of how to find ui as a function of the random jump times and 
the initial condition. Below it will be convenient to use to = and uq = 0; with that convention Eq. ( [Tl| ) 
holds also for i = Q if we take the special initial condition = X + a. 



2.2 Basic integral 

The persistence probability Q{T) can be expressed as a path integral on all random functions £,{t), hence 
as an integral on all jump times ti,t2, - ■ ■■ It is now useful to note that the ui are ordered according to 
= Mo < < U2 < . . . , so that there exists an L > for which 

UL<T < UL+i (12) 

The interpretation is that after the ith jump the function ^(t) is sure to stay above the threshold X, even 
if no further jumps occur, in the interval [0,T]. Summing on all possibilities implied by Eq. ( p^ ) we can 
write Q{T) as 



i^O JtL^l 



ul)Q{ul+i-T) (13) 



where Q is the Heaviside step function and where we used that e"''*^ is the joint probability density 
for the first L jumps to occur at ii, ^2, • • • , ^l- The L = term in Eq. ( p^ ) has no integrals and is equal to 
Q{ui — T). In the remainder we will use the shorthand notation 

//•111 rn-2 rut 

= / / dh dt2... dte £=1,2,... (14) 
Jo Jti Jtt^i 

The expression ( [l^ ) for Q{T) bears great similarity to the grand-canonical partition function of an assembly 
of interacting particles in a one-dimensional volume T, with the jump times ^1,^2, • • ■ playing the role of 
the particle positions and with the interaction implicit in the upper integration limits Mi,it2, . . . . 
In terms of Laplace transforms Eq. (|l^) is equivalent to 

/>OC 

Q{cu) = / dre-"^Q(T) 
Jo 

00 „ 

= uj-^Y. / e-''*^(e-'^"^ -e-""^+0 (15) 



L=0 " ^ 

One more rewriting is useful. For L > 1 one easily finds the relation 

L~l JL-1 

When Eq. ( |l6| ) is substituted in Eq. ( |l5| ) cancellations occur. After we replace cu with the dimensionless 
variable 

n = P±^ (17) 
p 
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we can express the problem by the three equations 

luQ{lu) = l-K{n) (18) 

OO 

K{n) = ^/^L(r!) (19) 
KLin) ^ J e-"''"^+i L = 0, 1,... (20) 

of which the last one implies, in particular, that KQ{n) = q-^p'^^ . Our task is to evaluate the phase space 
integral in Eq. ( pO[ ) and to find the relevant nonanalyticity of Q(w). In terms of the Laplace variable uj 
the persistence exponent 9 will be given by 

-uji^ -p{fii - 1) (21) 

where uJi is the real part of the rightmost nonanalyticity of Q{uj) in the complex co plane, and Qi is the 
corresponding value of O. Any further nonanalyticities at ^,2,^3, ■ ■ ■ will similarly give rise to correction 
terms in the decay of Q{T) characterized hy 62,03, ■■■ ■ 

2.3 Interaction potential V{y) 

At this stage the problem is to calculate of Eq. (|20|), defined as an integral via Eq. (p^), in which the 
upper integration limits ue are defined recursively via Eq. (p^). This problem depends parametrically on 
the function A{£^) or, equivalently, on f{t), and on the threshold X. We can still gain by transforming to 
another set of parameters. That will be the purpose of this subsection. 

Each jump provides the process with an additional lapse of time before hitting the threshold. The 
extra time furnished by the ^th jump is m^+i — Ui. The negative slope of / restricts u^+i — to a maximum 
value that we shall call r and which occurs for tg = u^. Using this in Eq. (|l|) we see that r is the solution 
of 

/(-r) - /(O) = a (22) 

where, of course, /(O) = X. The £th jump will generally take place before rather than at the ultimate 
instant ug. Due to the upward curvature of / the actual extra time gained is therefore generally less than 
T. We will express this curvature effect explicitly in terms of a variable vg by setting, for £ = 1,2,..., 

ug+i - ue = T ~ vi (23) 

whence necessarily < vg < t. We now use this equation in (O) to eliminate ug+i and we then subtract 



Eq. (|22|). This gives 

fiU -ue-T + vi)- f{U - ui) ~ f{-T) + /(O) = (24) 

from which Vi can be solved in terms of ui — tg. Although the jump density p does not appear in the 
above equation, it will turn out to be convenient to write the solution vi in the scaled form 

pvi = V{piut.~'ti)) (25) 

in which V has the expansion 

00 

V{y) = rY,gky' (26) 

fc=l 

It is easily seen that in accordance with Eq. ( p4[ ) one has vg — when ug — ti = 0. One obtains from 



Eq. (24) an equation for gu in terms oi gi, . . . , gu-i by differentiating k times with respect to ug — tg and 



setting ug — tg — vg = 0. This yields for the first two coefficients 



^ r(-r)-/(0) ^ r(0) _ /(0)V"(-r) 

r/'(-r) ' 2pr/'(-T) prfi-rf ^ > 
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We emphasize that we do not suppose r small. In cases where the limit r ^ may be taken, obvious 
simplifications occur. 

We continue now the analysis of the integral (^0|) for ■ This analysis may be performed for general 
initial condition ^o; however, from here on we shall impose ~ i, whence ui — r, in order to have simpler 
expressions, knowing that the persistence exponent will not depend on ^q. We will briefly come back to 
this point after Eq. It is useful to define r = pr, which, by Eq. (^2|) and relation (|l0| ) between A and 
/, is equivalent to Eq. (|3|) of the Introduction. Rewriting Eq. (^o|) in terms of the new integration variables 
ye = p{u£ — ti) and using Eq. (|2^) iteratively to express ul+i in terms of the yi we find, for L = 1,2,..., 



rW{yo) rW(yi) 

KL{n) = e-(^+i)"'- / dyi e"^(^i) / dy2 e"^(^^) 
Jo Jo 

>-W{yL-i) 

dyL e"^(^^) (28) 







where we have abbreviated 



W{y)^r + y-V{y) (29) 

and, by convention, put yo = 0. A special case is Ko{i^) = e~^^ . We have now transformed the phase 
space integral for to a problem depending on the parameter r and the interaction potential V{y). The 
original parameters X, a, and the function A(^) [or, equivalently, f{t)] no longer appear. 



3 Noninteracting theory: V = 

The noninteracting case is obtained by setting = in the preceding development. Strictly mathemat- 
ically it is not needed to study this case before passing to the next sections. However, from a physical 
point of view it is highly desirable to have a good idea of the noninteracting system before introducing 
interaction. 

For y = the theory depends on the single parameter r. Correspondingly, all derivatives of f{t) 
beyond the first one vanish and f{t) is given by 

fit) = X- f{0)t (30) 

We shall denote quantities referring to this linear decay curve by an index 0. When combining the above 
expression for f{t) with Eq. ( |2^ ) and the definition r = pr wc find that in this noninteracting case r is 
given by 

pa ,oT\ 

which is an instance of Eq. (^ with go — 0. 

Combinatorial problem. Having thus found the parameters of the noninteracting problem, we have to 
substitute them in the general expression ( ^8|) for KL{n). Imposing as before the initial value — X + a 
we obtain, after changing to the integration variables xg — yt/r^, 

pl+xo rl+xi rl+XL~l 

ifi(r!) =e-(^+i)^'-«rM dxi / dx2.../ dxL (32) 

JO Jo Jo 

where xq = 0. The L-fold integral in the above equation, that we shall refer to as II, constitutes the 
heart of the problem. In terms of the analogy with an L particle system the xe are the particle positions. 
There is no energy associated with the allowed configurations (xi, . . . ,xl), and log//^ is the entropy of 
the system. 

Upon converting to the integration variables se = £ — xi, where £ = 1, . . . , L, we have 

/■ 1 /'2 pL 

Il= dsi / ds2 . . . / dsL (33) 
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The same integral but with all upper integration limits set equal to L + 1 is elementary and equals {L + 
1)^/L\ It represents the phase space volume for putting L points on (0, L+l], not counting permutations 
as distinct. Hence II = [{L + 1)^ /L\]pL, where pl is the probability that L randomly chosen points on 
(0, L+l] are such that, for fc = 1, . . . , L, the number Mk of points in the interval {L + l — k, L+l] is less 
than k. 

This may still be rephrased as the following nonelementary combinatorial problem. Let L balls be put 
randomly in L + 1 numbered vases; then pL is the probability that the first k vases contain together at 
least k balls, for A: = 1, 2, . . . , L. 

We found no direct way to calculate pl and invoke a theorem due to Takacs, of which we adapt the 
proof to the present context in Appendix A. The result is that II — {L + 

Persistence exponent. Using this in Eq. (^2|) and substituting in Eq. (fol) we have 



(L + 1)^-1 

e J'o ■ 

L=0 



E 

L=0 



L\ 

g-L(aro-logro-l)+C'(logL) j^g^-j 



where in the last step we have used Stirling's formula. It is clear that as Q. is lowered, K{^) diverges when 
r2 attains a value that we shall call llo and which is given by 

f^o(ro)---log- + - (35) 
ro ro ro 

Because of Eqs. (^) and (|2l|) the persistence exponent is 

0o = p(-log- - - + l) (36) 
^ ro ro ro ' 

Converted to the original variables of the problem this becomes Eq. (^) of the Introduction. This exponent 
may also be arrived at in ways independent of the recursion relation formalism of this work {e.g. with the 
aid of the method of Ref. Appendix A), and appears in other contexts as well [e.g. the recent work 
of Bauer et al. ]l6t). It will appear again in the next section at the end of a very different calculation. 



4 Interacting theory: V > 
4.1 Small curvature limit 



The interacting theory has > in Eqs. (^8|) and (^9|). We will not be able to treat the general case, but 
only the one in which the series ( p6| ) for V{y) is dominated by its linear term. Curiously enough, although 
we have to suppose V small and although our final results for the exponent 9 will be perturbatively close 



to the zeroth order expression (36) of the previous section, the solution method of the present section is 
nonperturhative in the sense that we do not start from the = solution, and that in the limit y — > 
the method of this section ceases to work. 



The linear term dominates the series ( |26| ) for V{y) in particular in the following two limits. 

(i) g at fixed r, with rgi = g and rgk — o{g) for fc = 2, 3, . . . ; 

(ii) g ^ and r — s- with a fixed ratio g/r = c. 

In both limits the curvature parameter g tends to zero, and we shall refer to them as instances of a 
small curvature limit. The developments of the next two subsections are common to both limits. We set 
g = Tgi, which, by Eq. (^) and relation (|l^) between A and /, is equivalent to Eq. (^) of the Introduction. 
Retaining only the linear term in V we get 

V{y)^gy, W{y) ^ r + {1 - g)y (37) 
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Hence we have a theory with two dimensionless parameters, r and g; for g = Q the noninteracting theory 
is recovered. 

In the devefopments that foUow the higher order terms, suppressed in Eg. (|37|) , may be taken into 
account perturbatively to show that their effect is negligible to leading order. Throughout the present 
section, the discussion will be only in terms of the interaction constants r and that we shall consider as 
independent parameters. In Section we will return to the original variables of the problem. 



4.2 Recursion for K{Vt) 

If Eq. is substituted in Eq. ( p8| ) , it becomes possible to carry the L integrals out recursively for 
arbitrary g, as we shall show now. It appears that one needs auxiliary functions -ftT^^"'' and /C'") with 
rt = 0, 1, 2, . . . . These are defined by 

oo 

(38) 

in which for L = 1,2,. 



^0 

r+{'i--g)yL-2 (•'■+(i-s)i/t-i 

dj/L-i e"f^^-i / e"[i-(i-f)" (39) 

Jo 

and where we have the special case 

kI^'\u) = exp[-r!r(l - ff)"] 

= En (40) 



When Eq. ( |37| ) is substituted in the functions K and Kl of the preceding section, one sees that K = K^^'i 

and Kl 
relation 



and Kl — K^l \ Upon carrying out in Eq. ( |39D the integral on we find straightforwardly the recursion 



- bn [4"+/^ - Kl-i] (41) 

where 

I g-Or(l-3)" 

^" " l-(l_g)«+i 



Eq. (41 ) is valid for L — 1,2,... and n = 0, 1, . . . , and must be supplemented with the boundary condition 
(|^) at i = 0. Substitution of Eqs. (|^) and (|l]) in Eq. (||) yields for the iC'"' the recursion relation 

= hn -K]+En (43) 

The existence of this recursion relation is the key to the success of the present method. We remark that 
for g — Q the coefhcients &„ are undefined and the recursion does not exist; hence this solution method is 
nonperturbative. 

If we apply ( ^ ) to _ftr = K^'^^ and iterate n times, the result is 

K{n) = {BoEo + BiEi + . . . + BnEn) 

-(B1 + B2 + ... + Bn+i)K{Q) + S„+iX("+i)(r!) (44) 

where we abbreviated 

n-l 

Bo = l, Bn=Y[bj 71=1,2,... (45) 
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We examine now 6„, Bn, and £'„ for n ^ oo. Eqs. (^0[), (|4^), and (45) show that in that hmit 

5oo = Br, ~ Hoo = 1 (46) 

where 



Hence for n ^ cxd we obtain from ( [43|) an equation for if with well-defined coefficients. Using (46) one 
readily finds the solution 

if(°o)(17)^^--^-^(")=^-^(") (48) 

1 - 6oo - 1 

For n ^ oo we now replace if in Eq. ( ^4[ ) by if (°°) found in Eq. (^). Upon solving for K{fl) we get 

^ ^ li,, i?oi^o + + ■ • ■ + B,.E,. + [i?„+i/(l - »-^)] 

^ ^ So + Bi + . . . + B„ + [B„+i/(l - r!-i)] ^ ' 

= |f7|>l (50) 

Z^n=0 

in which i?„ is given by Eq. (Eol) and where from Eqs. (E5|) and (B5) we have i3o = 1 and 



Bn = ^ exp 



f7r 



1 - (1 - 5)" 



ni^TT^ n.1,2,... (51) 



J f \ 1 - (1 - g) 



Expression (|5C|) constitutes the solution of the problem of this work; the remaining analysis is needed to 
extract the persistence exponent 9 from it. Eq. ( |50| ) holds for the initial condition = X + a] without 
giving the proof we state that for general the same expression (|o|) is obtained except that in the 
definition (HO) of the En one should replace r by pui and remember that /(— ui) = ^o- 

By Eq. (pTl) we have 9 = — p(rii — 1), where fii is the rightmost nonanalyticity of K[fl) in the complex 
n plane. We expect the relevant nonanalyticities to be due to zeros of the denominator of Eq. (pO[), for 
which we shall introduce the special notation 



In view of the remarks of the preceding paragraph this denominator is independent of the initial condition 
^0- Obviously its zeros can occur only for 17 < 0. It is furthermore clear in advance that for g > 
the persistence probability must decay at least as fast as for g ~ 0, whence 9{r,g) > 9(r,0) — 9Q{r). 
Consequently we expect that fli < no{r), where Hq is the function defined in Eq. (pq). 

4.3 Analysis of K{n) 

We shall evaluate H(fl; r, g) asymptotically in the two limits — s- at fixed r, and 5, r — s- at fixed g/r. 
In order to prepare for these limits we will transform the sum on n in Eq. (^2|) into a contour integral, to 
which we shall then apply the stationary phase method. 

It is first of all necessary to extend the definition of the summand i?„ to arbitrary complex n. To that 
end we consider the function 

n,°:i(i-(i-gp ^ 
n^i(i-(i-5)^-^+^) 

which on the positive integers reduces to 



(^) = (53) 



^g{n) = X{{l~{l-9y) n = 2,3,... (54) 
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and Tg{l) = 1. This function Tg{z) was introduced in 1847 by Heine (see Ref. ||T^); nowadays it is 
usuaUy defined ||l^ with an extra factor g^~^ on the RHS of Eq. (p3|), and then caUed the q-gamma 
function^ where q — 1 — g. The function Tg[z) of Eq. ( |5^ ) has various properties reminiscent of the 
ordinary gamma function. In particular, it has poles for z = 0, —1, —2, . . . , and the residue Rm m z — —to 
is equal to 

fl„ = (-ir.g-i(l-ff)W'"'^) / m = 0,l,2,... (55) 

rg(TO + 1) 

We can now express H of Eq. ( |5^) with the Bn of Eq. (|5l| ) as 

Hin-.T.g)^^ [ dze^(^'") (56) 
27ri Jc 

in which 

h{z,n) = zlog(-r!)-f]rg-i[l- (1-5)^1 

^^z{z+l)log{l^g) + logTgiz) (57) 

and where the integration path encloses the poles of Tg(z). Equivalently, we may let this path run from 
—oo to below the real axis, encircle the origin counterclockwise, and run from back to ~oo above the 
real axis. The poles inside this contour exactly generate the terms of the series in Eq. (|5^). A factor (— 1)" 
coming from (— i^)" cancels against the (—1)" from 



4.4 Limit g ^ at fixed r 

If one now scales with g according to v = gz and writes 

h{z,n) ^-h{i^,n;g) (58) 
9 

then the limiting function limg^o h^v, 51; g) = hiy, £7) exists and is equal to 

h{v,n) = v\og[-^l)~nr{l~e'') + \v'^ + [ d^ log(l - e"'') (59) 

2 Jo 

The poles having become dense, this function has a branch cut along the negative real axis in the complex 
V plane. 

Stationary points. In the limit g — *■ we may apply the stationary phase method. It appears that 
h{v, f2) has two stationary points u±(yi). There is a critical value Sic such that for O > f2c the u± are real 
and positive, and for f2 < ilc they are complex conjugate. At Sic we have v- = = Vc- The values Sic 
and Vc are the solution of 

K(yc, Sic) = 0, K,{vc, Sic) - (60) 
where the indices on h indicate derivatives. These solutions are easily found and read 

Sic = log - + - 

r r r 

= + (61) 

We see that Slc(r) = Slo(r), which establishes the relation of this nonperturbative calculation with the 
solution of the noninteracting theory given in Section ^. The analysis can be refined in the vicinity of 
[vci Sic)- Upon performing a double Taylor expansion in 

Ai^ = i/ - i^c, ASl = SI - Sic (62) 
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we find, in obvious notation, 

h{i^, n) = h{i^c, ^c) + Anhn + ^Ail'^ hnn + Ai^Ail h^n + ^ Ai^^ + . . . (63) 

wliere all derivatives are evaluated at (fc, ilc), we have used that h,^ = huu = 0, and the dots indicate the 
remaining third and the higher order terms. The derivatives that it will be useful to know explicitly are 

hun=r, h^^^ = -\og^r, /in = -r (log(l + ^^)) + logr) (64) 

The stationary point condition dh/dAv = applied to Eq. ( |63| ) now shows that Ai/ has to scale as Ail^/^ 
and we find 

i^±(n)^i'cT 7^ (2r Arj) (-g5^ 
logr 

The stationary point integrations are easily carried out. For AQ > the relevant stationary point is i/_ 
and the outcome of the integration is positive. For ACl < the complex conjugate points both contribute 
and the result is 

H{n; r, g) ^ 2Ho{n; r, g) cos (- f^^A^A ^'^ + -\ (66) 

\g\9h,y^„ I 4/ 

where HQ(fl;r^ g) is positive and where it should be remembered that and Afl are negative. 

Zeros of H(ri;r, g). Upon substituting in Eq. ( |66| ) the explicit expressions Eq. ( p^ for the derivatives 

of h we see that the function H has zeros for ft = ilj with 

1 / l\2/3/ \2/3 

^i=^c-^[\og-) ((4j-l)— g) J = 1,2,... (67) 

For j ~ 1 we obtain the rightmost singularity of K(Q) in the complex Q plane. Hence by Eq. ( ^ ) we 
obtain for the persistence exponent 9 the result 

1 / l\2/3/q7r \2/3 

9 ^ Ooir) + p-[\og-j i^—gj r fixed, 5 ^ (68) 

with the function 9o given by Eq. (|36|). When reconverted to the original variables this gives the result 
announced in the Introduction. The second term on the RHS of Eq. ( |6^ ) represents the leading order 
curvature correction to the persistence exponent. It is nonanalytic at zero curvature. 



4.5 Limit g,r ^ with fixed ratio g/r 

We set g = cr. This limit requires an independent evaluation starting from Eqs. ( p6| ) and (p7|). The critical 
point (j^c, r^c) is still given by Eqs. (61), but it appears necessary now to scale the deviations from it as 

AD = log^ - Av, An^r log2 - An (69) 
r r 



Setting for convenience = 1/ log(l/r) and expanding in one finds, instead of Eq. (|63|), the expression 
h{iy, n) = h{iy„ n^) + (el + 2et)An + e^AUAD - ^4AD^ + O(e^) (70) 

in which /i(fc,ric) is itself of order Cr. The stationary point condition dh/dAv now leads to AD± — 
=F(2Afi)i/2, which when the original scaling is restored is the same as Eq. (|65|). The integration through 
the stationary point involves only the two 0{ef) terms in Eq. (|70|). In view of the proportionality between 
g and r the curvature in the stationary point is in this case of order {glog'^ g)~^ instead of g~^. After the 
calculation is done the expression for 9 appears to be exactly what one obtains by naively substituting 
g = cr in Eq. (|68|), that is, 

o /97r \ 2/3 1 / 1\ 2/3 
+ -l,(logi) (71) 



Eqs. (|6g) and ( [7 If ) constitute the main result of this section. In the following section we shall compare them 
to direct Monte Carlo simulations of the process S,{t). The zeros ^2,^3, ... , whose explicit expression is 
furnished by Eq. (|67|), lead to exponentially small additive corrections to the leading decay of Q(T). 
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5 Examples 



In the foUowingapplications wc will start from functions A{^) defining specific examples of the Langevin- 
type equation (|l|). 

First example. If in Eq. (0) we take A{£_) = the result is the linear equation 



dm 

dt 



(72) 



The parameters r and g follow directly from Eqs. (^) and (|^), respectively, with the result 



f log(l + ^) 



9 = 



X 



As X becomes large, r and g tend to zero simultaneously with the fixed limiting ratio 



.9 /3 
c = lim — = — 

x^co r p 



(73) 



(74) 



Hence we are in the situation of Subsection 4^. Eq. (22) then leads to r = (3~^ log(l + a/X), so that for 
X — > oo we have r — > 0; since y in Eq. (Um is of order r, the series for V{y) is one in ascending powers 
of T and we were justified in Subsection |4.l| to neglect the nonlinear terms in V{y). If we now substitute 
expression (fz^ ) for r and c in Eq. (Iti] ) and neglect subleading terms in the curvature correction, we find 
the level exponent 9{X) given in Eq. (||) of the Introduction. 

We have performed a Monte Carlo simulation of the proces ^(<) of Eq. (Q and determined the per- 
sistence exponent 0{X) for X ranging from the average {£,) = pa/ (3 up to eight times that value. Fig. 1 
shows the Monte Carlo data for 9{X) along with the theoretical result, Eq. (^), for asymptotically large X. 
There are no adjustable parameters. The dashed curve ("free theory") represents only the first term on 
the RHS of Eq. (||); the sohd curve ("interacting theory", full Eq. (||)) includes the leading order curvature 
correction, which is the main result of this work. This correction appears to be an important effect. The 
excellent agreement between the interacting theory and the simulation data strongly suggests that higher 
order corrections to Eq. (0) vanish for X — oo. 



200 




Figure 1: Persistence exponent 9 as. a. function of the threshold X for the process ^(t) of Eq. (|7^ ) with 
a = 1, p = 10, and (3 = 2. The average value of this process is (^) = 5. The error bars of the simulation 
data are smaller than the symbols. The interacting theory expanded to leading order (solid line, Eq. (||) 
of this work) is in excellent agreement with the simulations. 
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Figure 2: Same as Fig. 1, zoomed on thresholds X close above the average (^) — 5. 



Fig. 2 shows a zoom on values the leading order behavior of the interacting theory (solid 

curve) still represents a considerable improvement over the free theory, but as X — > (^) , higher orders in 
the expansion become necessary. For X < (^) the expansion of this work does not apply. 

Second example. Let A{^) be such that for some small parameter e 

A{0 = A{eO (75) 



Eq. (10) may then be recast in the form 



^t = - ^ (76) 



X 



Aix) 



whence it follows that f{t) scales as 

f{t) = e-'T{et;eX) (77) 
We have by construction /(O) = X as before. Furthermore 

/(")(0)^efc-i^(fe)(0;eX) (78) 

where the differentiations of J- are with respect to its first argument. If now we agree to choose X of order 
e^^, then the kth derivative of / is of order e*"'"^. This guarantees that gk is of order e'^, as is illustrated 
by Eqs. ( p7| ) for gi and g2- Hence the conditions of limit (i) are fulfilled and the calculation of Subsection 
4.4 apphes. 

In order to test the nonanalytic dependence on the curvature parameter g in Eq. ( |68| ) we have performed 
a Monte Carlo simulation of ^(t) for the particular choice 

f{t) ^X~ I3t + ef (79) 



that is, T{x, y) — x ~ f3y + y . The corresponding A follows from Eq. ( [76| ) and by means of Eqs. (|75| ) and 
(^) we find 

/ 4fr7\-i/2 

5 = l-(l + ^) (80) 

Monte Carlo simulations were carried out at fixed a and (i for various values of e. In Fig. 3 we show the 
persistence exponent 6* as a function of 5, together with the theoretical g^/^ law of Eq. (|68|). The agreement 
is excellent. 
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Figure 3: Persistence exponent 6 for fixed threshold X as a function of the curvature parameter g for the 
second example of Section]^ with a = 1, p = 10, and j3 = 20. The solid line represents the theoretical g^/^ 
law (Eq. (|6^) of this work) , which appears to provide an excellent description of the data. 



6 Limit of Gaussian noise 



The Langevin equation (with white Gaussian noise) and its extension to colored Gaussian noise are at 
the basis of much recent work on persistence; see e.g. the recent review by Majumdar Q. There is a 
large body of knowledge today about the persistence properties of such Gaussian Markovian processes, 
and a perturbative method around the Markovian case has recently been devised by Majumdar and Sire 
[Q (see also Oerding et al. and Majumdar et al. Q). The equation of this work, Eq. (|l]), with jumps 
of arbitrary finite size a, provides, on the contrary, an example of strongly non-Gaussian noise. In this 
section we show how for a — )■ the Gaussian limit is approached. This limit, just as the one of zero 
curvature considered in Section ^, is a singular point in parameter space. 



6.1 Gaussian persistence exponent 6g 

Let C,(t) obey the linear Langevin equation 

^=-/3C + i(0 (81) 
where L{t) is Gaussian white noise of average {L{t)) = and correlation 

{at)at')) = TSit-t') (82) 

The level exponent 9g{Z) for this process, associated with the probability for C(t) not to have crossed a 
preestablished threshold C = Z in a time interval has not to our knowledge been calculated in the literature. 
The related exponent associated with crossing upward through the threshold has been considered by 
Krapivsky and Redner (see also Turban jl^). It is easy to find Og{Z) by a method similar to theirs, 
as we will show now. The probability distribution P(C, t) for the process ( |8l] ) evolves in time |l^ according 
to the Fokkcr-Planck equation 

dP d F 

The persistence exponent Oq is the eigenvalue of the slowest decaying mode for Z < C < oo satisfing the 
boundary condition P{Z, t) = 0. We set P{C, t) = -P(C) exp(— 0g t). It is well-known that the equation for 
P{C) can be transformed to the eigenvalue equation for the quantum harmonic oscillator. This fact has 
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been exploited in previous work ^, ^ on persistence exponents. Here, in view of our interest in the 
interval Z < < oo and the limit of large Z , we must transform 

P(C)=c^^V(C) (84) 
with C = X^/^{y/2P/rC - 2VA) and a = 1/2 + 6'g//3- Then P(C) satisfies the eigenvalue problem 

^ - C (l + ) P(C) = 0, P{Z) = (85) 

In the limit of high threshold Z we expect 9g, and therefore A, to diverge. Hence in this limit 

^ - C) PiC) = 0, P{Z) = (86) 

The solution of Eq. (^) that vanishes for C — > cx3 is the Airy function Ai(C). The boundary condition 
Ai{Z) = leads to Z = ai, where ai = —2.3381.. . is the first zero of Ai. This condition fixes Oq in 
terms of Z; upon expanding for large Z one finds 

eG(^)-^ + |ai|/3i ^ +... (Z^oo) (87) 



2r ' V 2r 

which is the desired result. 



6.2 Gaussian limit 
6.2.1 Limiting procedure 

In Eq. (|l]) we substitute now ^ = C + paj (5 and a^j, 6{t — tk) — pa ^ L{t) and take the "Gaussian" limit, 
defined as 

p — !■ cx), a — !■ with r = pa^ fixed (88) 

The result is that Eq. (^l|) appears. One easily verifies that {L{t)) = and that the cumulants of L, which 
for n = 2, 3, . . . are given by 

n-l 

{L{h) . . . L{tn))c = pa" n ^^^^ - (89) 

k=l 

vanish in the limit of Eq. (|8^) when n > 3. Hence L{t) is Gaussian white noise. The above transformation 
changes the threshold X into Z = X — pa/ (3. One now expects that the Gaussian persistence exponent 
9g{Z), found by direct calculation at the end of the previous subsection, should also be accessible as a 
limiting case of our general approach. Naively, one may attempt to obtain 9q (Z) by taking the Gaussian 
limit, followed by the limit Z —> 00, in expression ( |68| ) for 0. After a short calculation that procedure 
leads to 

This differs from the exact result, Eq. (|87|), only by the numerical value of the coefficient of the subleading 
term; moreover, the difference ((97r/8)^ ~ 2.3203 . . . versus |ai| — 2.3381 . . .) is only about one percent! 
Nevertheless, Eq. (^) is right and (^0|) is not. The rather obvious reason is that the Gaussian limit (which 
implies aZ — > 0), followed by Z — s- 00, does not commute with the limit that was taken to arrive at 



Eq. (68) {viz. r, 5 — > at fixed r/g, which implies aZ — > 00). In order to find Oq{Z) within the formalism 



of the preceding sections it is necessary to start again from the integral representation of H(p,\r,g) in 



Eq. (56). Below we will see how to do that. 
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6.2.2 Calculation of 60 

Let us consider H{fl;r,g) of Eq. js^). In view of Eqs. (p7|)-(p9|) it is represented as an integral on v of 
the function exp((/~^/i(j/, The Gaussian limit is controlled by the parameter a, which should tend to 
zero. At fixed T = pa'^ and Z = X — pa/ (i we find from Eq. ( [73| ) that in that limit g = 70^ + O{o?) and 
r = 1 - + 0{a?) with 7 /3/r. We recall now Eq. (|2|), which says that 6 = p{l - Expecting 9 
to approach a finite limit Og, we set J7 = 1 — a?W ^ where W is the appropriately scaled variable for the 
relevant region of the complex frequency plane. Hence, if the rightmost zero of H(Q]r,g) in this plane 
occurs for W — Wi, then 

0G = lim p{l - rii) = rW^i (91) 

Stationary points. As a preliminary we consider the stationary points of h{iy,n). Expanding the 
equation dh/dv = for small a while anticipating that e" will be small we find that these points are 
solutions of 

- a^W + (1 - a-iZ)e'' - - ^e^"" + . . . = (92) 

where the dots represent terms of higher order in a and e"^. This shows that there exist solutions with the 
scaling Re ^ logo for a ^ 0. Solving explicitly we obtain 




a^Z\~l±Jl-—-.\+Oia') (93) 



In the above expression there appears a critical value of W equal to Wc = ^■y^Z^. For W>Wc, which we 
expect to be the relevant regime, the stationary points therefore are i'± = —A — in ± ip* with 

-A = log{aV2W) + 0{a log a), p* ^ arccos -^L + 0{a) (94) 

Instead of the variable of integration v we will henceforth use p defined by 

1^ = -A -m + ip (95) 
We will not exploit directly, in what follows, our knowledge of p* . 



Gaussian limit. We consider h{i', fl) of Eq. (59) as a function of p. After some calculation we find that 
for small a 

h{v, n) = h{^A - iTT, Vt) + a^-fkip, W) + 0{a^) (96) 

with 

ryk{p, W) = ~\Wp + 7ZV2W^(e^^ - 1) - ivK(e2^^ - 1) (97) 
In the limit a — > the fimction H{Q;r, g) may therefore be rewritten as the integral 

H{n; r,g) = D J dp e'^^"''^) (98) 



with k(p, W) given by Eq. ( |97| ) and where D diverges when a goes to zero. However, D will divide out in 
Eq. ( [50| ) against the same factor in the numerator of K(Q). This completes the Gaussian limit. There is 
no small parameter left in the integral in Eq. (p8|). 

Limit of large Z. This integral may be reduced to a more elementary one in the limit of large threshold 
Z. The reason is that then the relevant values of W are close to Wc- We adopt the scaling 

W = PFc(l + = ^j^Z^{l + wZ'^/^) (99) 
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which wiU be justified by the results. We now consider the fuU Taylor series in /i of k{iJ,,W). Upon 
expanding each of its coefficients for large Z and retaining only the leading term we get 



k{fi,W) = -IjnjZ^/H^i-l^wZ^/'^'^'^^ 



2 ' 2 ' 2! 

n— 3 



-7Z2 ^(2"-2_i)Wl^ (100) 

n— 3 

If now the integration variable is scaled according to /i = A(7Z^)~^/^, then in the large Z limit all terms 



in Eq. (|100| ) except those with n = 1 and n = 3 go to zero. We are left with 

F(17;r,g) - /dA e-^T'^'^i^+^i^' (a = 0; Z cjo) (101) 



which is the integral representation of the Airy function. The only dependence left is on the variable w. 
Let the rightmost zero of H(fl;r,g) in the complex frequency plane occur for w — wi. We see now that 
wi is the solution of Ai(7^/'^wi) = 0, whence 

= |ai|(2r//3)2/3 (102) 

Upon relating wi to Wi by Eq. ( |9^ ) and using Eq. ( |9l] ) we finally get the expression of Eq. (|7|) for Oq- 

Discussion. It is instructive to return to the quantity fi* given by Eq. (^). The two stationary points 
are separated by a distance 2/z*, and substituting the various scaling transformations we see that, as 
Z oo, they have in terms of A the finite distance 2A* = 27-*^/^^-'^/^. We now observe the mechanism 
that is at work here. In Section ^, for a finite, hence far from the Gaussian limit, H{p.;r^g) is the sum 
of contributions from two stationary points at infinite separation (~ with 3 — > 0) in the v plane; as 
the Gaussian limit is approached, the two stationary points come within finite distance of one another, 
and their contributions cannot be separated any longer. This " interaction" between the stationary points 



leads to the replacement of the cosine in Eq. (B6[) by the Airy function in Eq. (101), and finally affects by 



about one percent the coefficient of the subleading term of the persistence exponent. 



7 Conclusion 

Beside many Gaussian persistence problems, there are also non-Gaussian ones occurring in statistical 
physics. We have pointed out and studied one class of such problems, associated with the specific non- 
Gaussian stochastic process that satisfies Eq. (||) . Its relation to several questions in statistical physics has 
been indicated in the Introduction. The sample functions of this process are deterministic curves inter- 
rupted at random instants of time by upward jumps. Among these, a zeroth order subclass is constituted 
by "random sawtooth" functions, characterized by linear decay with fixed slope. The persistence exponent 
Oq of this subclass is easy to find. We then perturb around this zeroth order problem by introducing in 
the decay a small curvature of strength controlled by a parameter g. As a consequence we have to deal 
with what is essentially a one-dimensional interacting particle system with coupling constant g, and the 
mathematics becomes considerably more complicated. The case of greatest importance covered by the 
present work is the linear equation, with exponential decay curve, that prevails for A{^) = fi^ in Eq. (^. 
Our result for this case is an asymptotic expansion, Eq. (^), of the persistence exponent 9[X) in the limit 
of high threshold X . 

The same equation for level X = {£,{t)), which is outside of the domain of the asymptotic expansion 
of this work, has recently been considered by Deloubriere . It would be of definite interest to extend 
Eq. (|^) to random upward jumps Ok at time tk, given that specific distributions of jump sizes Ok naturally 
occur in several models of statistical physics [ill . 
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A A theorem by Takacs 

We consider the problem of determining the probability pl that occurs Section |^. Let the variables Mk 
be those defined there. It is natural to set in addition Mq = and M^+i — L, so that our problem is to 
find 

PL = Prob{fc~ Affe > for fc = 1,...,L+1} (103) 

Relevant to this problem is Theorem 3 by Takacs [|l9| , which concerns nondecreasing random functions 
on line segments. The author jl^ indicates that this theorem has an analog valid for nondecreasing random 
sequences. For the present case the full proof runs as follows. 

The range of the index k may be extended to arbitrary positive k by the definition 

ML+i+k = L + Mk (104) 

This amounts to repeating the set of random point onO<s<L + l periodically in the segments 
n{L + 1) < s < {n + 1){L + 1), where n = 1,2, ... . The random variable M^+e — Mk, where fc = 0, 1, 2, . . . 
and £ = 1,2,..., represents the number of points in the interval k < s < k + i, and the probability 
distribution of this variable is obviously independent of k. Let now for fc = 0, 1, 2, . . . 

r 1 if Mk+e-Mk <^ for £= 1,2,... 
4 = <^ (105) 
[ otherwise 

Then the probability distribution of Sk does not depend on k, and Sk+L+i — ^k- It is easy to verify that 



Mk+i - Mk < i holds for all £ if it holds for £ = 1, 2, . . . , L + 1. Hence Eq. ( |103D shows that pl is the 
probability that So be equal to 1. We may write equivalently pl = (<5o), where the average is on all random 
sequences Mi, . . . , M^. But since all 6k have the same distribution, hence the same average, we also have 



P^ = T^l = (E Sk) (106) 



k=l k=l 

We consider now the sum on the 5k in the last member of the above equation. The condition for 6k to 
equal 1 may be rewritten as 

j-Mj>k-Mk for aU j = fc + 1, . . . ,fc + L + 1 (107) 

In the range k<j<k + L + l the function j — Mj has the initial value k — Mk and the final value 
k + L + 1- Mk+L+i ^ k- Mk + l, where we used the definition (|104[) . If 4 = 1 , then j - Mj > k~ Mk + 1 



for all J = k + 1, . . . , k + L + 1, and this means that 6k+i — . . . — 6k+L — 0. Hence '^^^ be 

equal only to or to 1. We now prove that in fact it equals unity. For it to be zero, all 6j in the range 
of summation would have to vanish, whence we would have 6j = for all j > k. There would then exist 
an increasing sequence {jr}^o (where jo = k) such that the corresponding sequence {jV — Mj^}^^^ is 
nonincreasing. This however is in contradiction with the fact that j — Mj increases by 1 whenever j is 



augmented by X + 1. It follows that X]j=fc ^k — 1, whence by Eq. (106) we obtain pL ~ l/(£ + 1) and 
Il = {L + 1)^-^/L\ 
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